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We investigate the behavior of flux hnes in the presence of both columnar and 
point disorder in 1+1 dimension using Renormalization Group (KG) and world-line 
Quantum Monte Carlo (QMC) techniques. In particular, we calculate the transverse 
wandering correlation function for a single boson and recover known results for point 
and columnar disorder separately. We then examine the existence of a localization 
transition of a flux line in the simultaneous presence of both types of disorder. We 
also performed RG for interacting flux lines in the presence of both disorder. RG 
indicates that the vortex glass is unstable with respect to arbitrary small amount of 
columnar disorder. Using QMC we find that the Bose glass transition temperature 
is reduced by the point disorder, in agreement with recent RG calculations. Further 
we find that the region posited to be an Anderson glass is completely destroyed. 
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I. INTRODUCTION 



There has been a great deal of interest lately concerning the creation of a deliberately 
disordered environment as a method to pin flux lines and create a true superconductor with 
a vanishing linear resistivity. Two types of disorder have been employed: point disorder 
which is uncorrelated in space, and columnar disorder produced by heavy ion irradiation 
which is correlated in a preferred direction [^. While much theoretical attention has been 
lavished on the behavior of flux lines in the presence of either point or columnar disorder 
and their respective pinned glassy phases, the interplay of the two types of disorder is still 
an open topic. Both types of disorder lead to pinning of flux lines at low temperatures, but 
the nature of the phases is quite different. In the vortex glass the flux line is pinned in a 
tortuous meandering path by the point disorder. In the Bose glass 0,^, the flux line remains 
relatively straight, locked to the columnar pins. However the interface of the Vortex and 
Bose glass phases, where point and columnar disorder compete, remains beyond the reach 
of present analytical approaches . The only guidance comes from a recent study, where it 
has been argued that a single flux line will always be localized for arbitrary weak columnar 
disorder, irrespective of the strength of the point disorder Interactions between lines 
could then further limit wandering and enhance the localizing tendencies. 

In this paper we analyze the localization of single and interacting flux lines in the pres- 
ence of both point and columnar disorder. In particular we employ a mapping of the flux 
line problem in d+1 dimensions onto the problem of disordered quantum bosons in d di- 
mensions. We measure transverse flux line wandering of the single flux line, and discuss the 
possible existence of a critical ratio of the columnar to point disorder strength to localize 
the boson. Turning then to the interacting boson problem, we study the phase boundary 
of the Bose glass phase and its stability to point disorder. We combine these simulations 
with a Renormalization Group (RG ) treatment of a disorder potential which interpolates 
between point and columnar disordered potentials and find evidence for trajectories scaling 
away from the vortex glass phase. We also find evidence that a localized phase at weak 
boson repulsion, the "Anderson-glass" region, is destroyed by the de-phasing effect of the 
point disorder. 



II. FORMALISM 
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A. Mapping 



Our Monte Carlo analysis is based on the mapping of flux lines in d+1 dimensions 
onto the quantum problem of interacting bosons in d dimensions. This mapping can be 
understood qualitatively as arising from the similarity between the flux lines aligned parallel 
to the applied fleld and the boson world lines which appear in a path integral representation. 
More specifically, the classical De Gennes - Matricon free energy in the absence of disorder 

Fciass= ^dz\J:h{dr,{z)/dzy + Wyi\r^{z)-r,{z) |)|. (1) 

Here the displacement field ri{z) describes the position of fiux line z in a direction perpen- 
dicular to that of the applied field {z). The interaction potential V between the fiux lines is 
often simplified to a delta function that prevents fiuxes from occupying the same site. Then 
one can define a discrete field that takes on only integer values 

A{x,z)=Y.Q{x~r^{z)), (2) 

n 

where G is the Heaviside step function, and rewrite the classical free energy as 

Fciass = J dx£dz f^^iiidA/dzy + j{dA/dxf^ . (3) 
The partition function in the canonical ensemble is then 

Z^iass = f D[A]e-^^^-/^. (4) 



Consider now the quantum mechanical free energy for bosons in d spatial dimensions. 
In the continuum limit 

Fquantum = cItJ dx i^- -{^^ d'^ / dx' + C.c) + - | (5) 

plus terms which can be absorbed into a chemical potential. \l/ is the second quantized field 
representing the bosons. Following Haldane [^], we can write the Bose density operator 

N{x) = ¥{x)^{x) = [d,B{x)] e-''"^("), (6) 

m 

where d^B^x) = Nq + ti[x) is a field chosen to ensure the discrete nature of the Bose field. 
Here A'^o is the average boson density and H is a small fiuctuation. Linearizing in H we 
obtain 
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Fguantum = dr J dx I —{dB/drf + -{dB/dxf , (7) 

where we discard nonlinear terms representing commensuration effects. Tlius we find that 
the partition function 

Zquantum = J D [Q] e-^-'"-™/^ (8) 

Comparison of the two partition functions leads to the analogy 
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and all the energy scales get divided by T. The same holds true when one includes the 
chemical and disorder potentials. 

B. Simulation Algorithm 

In this paper we consider a discrete lattice version of the interacting boson problem, the 
Bose-Hubbard Hamiltonian |^ 

H = -tj^i^+io-i + + +^niec{i)- (9) 

i i i 

As we shall see, ec represents columnar disorder. The inclusion of point disorder will be 
discussed below. 

The partition function Z can be written by discretizing the imaginary time (3 = LAt, 
separating the exponentials of the kinetic and potential energies, and inserting compete sets 
of occupation number states: 
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Z = Tt e-^^ = Tr [e-^^^e"^^^]^ 

= E (n(^,l)|e-^-^e-^-^|n(^,2))(n(^,2)|e-^-^e-^-^|n(^,3)) 

n(j,r) 

. . .(n(z,L)|e-^"^e-^"^|n(z,l)) (10) 

The sum is over a classical occupation number field whose first index runs over the spatial 
sites in the lattice and whose second index runs from 1 to L and labels the "imaginary time 
slice", that is, the position of the complete set in the expression for Z. Since H conserves 
particle number, and K moves bosons only locally, the allowed configurations of n{i, r) form 
continuous world-line trajectories which are the analogs of the flux lines. The sum over all 
configurations necessary to evaluate Z is performed stochastically. 
Since P is diagonal in an occupation number representation, 

Z= J2 e"^' (^(^,l)|e-^"'^|n(^,2))(n(z,2)|e-^"^|n(2,3)) 

n(j,r) 

. . .(n(z,L)|e-'^-^|n(z,l)) 
Sb = VAt n{z, rf + Ar ^ e,{i)n{z, r) (11) 

i,T i,T 

The matrix elements of the kinetic energy operator K can be evaluated using the checker- 



board decomposition. |T0| We are now in a position to see how "point disorder" can naturally 
be incorporated as an additional imaginary time dependent random term in Sb which takes 
the form 

S',=ArJ2^pihrMt,T). (12) 

All the standard methods of world line simulations can now be employed. In addition 
to Monte Carlo moves which locally distort the world lines, we also include moves which 
shift the entire world line position spatially. This significantly helps equilibration in the 
presence of columnar disorder. The point and columnar disorder are bounded and uniformly 
distributed, i.e., —c,—p < ec,p < c,p, where c,p stands for columnar and point disorder 
strengths, respectively. 

III. SINGLE FLUX LINE 

In this section we consider a single boson in both a random point and columnar disorder 
potential and investigate the localization of the boson as a function of the strength of the 
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respective disorders, c/p. In particular, we focus on the transverse wandering correlation 
function. 



i?(r) = ((x(r)-x(0)f), (13) 

which measures the transverse extension of the boson world line from its initial position as 
it traverses the imaginary time axis. It is well known that in the absence of disorder the 
boson dynamics is diffusive, with the diffusion constant governed by t, or equivalently, the 



flux line tension [|1T| such that for large r, 

i?(r) = t T, no disorder. (14) 

Point disorder has been shown to promote flux line wandering due to the presence of favorable 
site energies which encourage excursions, and consequently, 

R{t) ~ r^'' point disorder, (15) 



with the wandering exponent C = 2/3 [|T2|. Lastly, in the case of columnar disorder simple 
quantum mechanics leads to the localization of the flux line for arbitrary pinning and thus 
the transverse extension of the flux line is bounded 

R{t) ~ Ic columnar disorder. (16) 

Since our boson must satisfy periodic boundary conditions, we study the wandering of the 
world line at imaginary time r = (3/2 after positioning the boson at its lowest energy site 
at time slice r = 0. Our runs use on average 100 realizations of the disorder. We first show 
our results which recover the above limiting behaviors in Figs. 1 and 2. For the case of no 
disorder (Fig. 1) we recover both the exponent 1 and the value of the diffusion constant t 
for the case of 32 site chain, in agreement with |TT|. The finite size of the chain limits the 



correlation function for large values of hopping t and for large imaginary time separation 
(3. For the case of point disorder (again for a 32 site chain, hopping t set equal to 1), we 
find that the wandering exponent C, increases from 1/2 to 2/3 for arbitrary point disorder 
strength, which is shown in Fig. 2 for p = 10 and 25. 

Turning on the point disorder p with a fixed value of columnar pinning strength c, we 
see in Figure 3 that R lifts up at higher (3 reflecting the nature of the point disorder which 
encourages flux line wandering. Our behavior is consistent with that of Kardar [O who 



via another method found numerical evidence, that for c/p > 0.83, the flux line is localized. 
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However, it has been shown that a simple extrapolation of the localization length to infinity 
is inadequate and the presence of rare regions due to the point disorder leads to slow glassy 
dynamics which keeps the localization length large but finite [Q. In Fig. 4 we show the value 
of the "localization length," i.e., the saturation of the wandering correlation function R, for 
runs at different values of c/p which indeed seems to indicate a crossover to delocalized 
behavior at a value consistent with that of Kardar. However, our simulations are limited 
to f3 values which cannot provide a closer analysis near {c/p)criticai to reveal whether the 
localization length slowly turns over and remains finite at lower values of c/p, as argued in 
Ref. i. 

IV. INTERACTING FLUX LINES 

We have seen in the previous section how the interplay between point and columnar 
disorder affects the dynamics of flux lines in the absence of interactions. We now turn on 
the flux line repulsion V to study the various pinned phases of interacting flux lines. The 
advantage of using a boson world line approach is that we can take the on-site interactions 
into account using the same numerical techniques as before. However, since the boson world 
lines are indistinguishable we found it more useful to measure the superfluid density ps 
rather than the transverse correlation function as an indicator of the Bose glass {ps = 0) 
and either the Vortex glass or Line liquid {ps 7^ 0) phases. In particular, we extend the 
results of our previous study of spatially disordered bosons by adding point disorder [jT4[. 



It was shown in Ref. 0] that there were two localized regions, separated by a superfluid 
phase as a function of boson repulsion. It was posited that the physics of the region at 
small repulsion (Anderson Glass) was qualitatively different than the phase at larger values 
of the repulsion (Bose Glass). Therefore the two phases may behave differently as point 
disorder is introduced. Localization of noninteracting bosons in ID is a consequence of the 
phase coherent (back) scattering of bosons and therefore the presence of a time-dependent 
potential will de-phase the boson wavefunctions and thus allow the boson to propagate |]15| . 



However, in the case of strongly interacting bosons, the boson repulsion governs localization 
and it could be that the effects of point disorder are suppressed. 

Recently, a renormalization group (RG) approach has been applied to the problem of 
flux lines with point and columnar disorder. This technique is able to capture the Bose and 
Vortex glass phases as well as the Line liquid phase which exists at high temperatures. It 
was shown that the location of the Bose glass transition is shifted towards higher repulsions 
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(i.e. lower temperatures) by the point disorder while the Vortex glass phase is unaffected by 
the columnar disorder. Also, a phase diagram has been constructed based on RG, suggesting 
a Vortex Glass to Bose Glass transition at a critical value of c/p. However, since the RG 
treatment was confined to small disorder, this transition could not be accessed in Ref. 0. 

We consider two approaches towards understanding the nature of the Bose and Vortex 
glasses in the presence of both point and columnar disorder. First we will reconsider the RG 
approach from Ref. |0 where we introduce a disorder potential which interpolates between 
columnar and point disorder potentials. This introduces another small parameter which 
allows us to construct the scaling trajectories. Similar ideas were used by Weinrib and 
Halperin on a related problem ||T^ . Then we perform Quantum Monte Carlo simulations to 
test the results of the RG analysis. 



A. Renormalization Group Approach 

We reconsider the coarse-grained free energy in a continuum description |]^, 

dxdz{ — {dxuf + —{dzuf — Vq{x, z)dxU — 2Vo(x, z)pcos{2n[xp + u{x, z)])}, (17) 

where u is the flux displacement-like field, p is the average line density, U is the on-site flux 
repulsion and ii is the line tension. The disorder potential Vq we take to represent both 
point and columnar potentials by writing 

{Vo{x, z)Vo{x', z')) ~ A6{x - x')h{\ z - z' \), (18) 

where we take h{z — z) = ^;qq^iyp, where zq is a short distance cutoff. The Fourier 
transform of h{z) is given by 

KK) ~ Co + cifc;(i-'^) (19) 

for small kz and co,i are positive constants. Thus we see that if a > 1 the potential is governed 
by short-range correlated disorder and thus describes point disorder, while for a < 1 the 
potential is governed by long-range correlations. Thus disorder with long-range correlations 
which decay at larger distances faster than 1/z leads to the same critical behavior as that 
due to short-range correlated disorder ||16|| . 

Replicating the above free energy and performing the disorder average leads to 

dxdz ^ djUa ■ djUj3 

— j dz'h{\ z — z' I) [A^9a;Ma(x, z)dxUp{x^ z + z') + cos(27r[Ma(x, z) — Up{x, z + z')])]^, (20) 
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where a, /3 G {1, ■ ■ ■ n} are replica indices and j = {x, z}. The bulk and tilt moduli are given 
by Kx = U/T and Kz = ci/T, respectively, and the bare coupling constant of the nonlinear 
cosine interactions is g = p^A/T^. Lastly, A^. = A/T^. 

We have performed a renormalization group analysis of Eq. (20) using 6 = 1 — a and 
g as small parameters. We developed a technique which is an extension of the sine-Gordon 
scahng theory as discussed e.g. by J. Kogut [|l^. We performed our calculations up to 
next-to-leading order in the couplings. This involved considering terms with 8 displacement 
fields. High scaling dimensions eliminate many of the newly generated terms. As a novelty, 
in this order the renormalization of the disorder correlation exponent a became necessary. 
Upon rescaling by a factor e', tedious but straightforward calculations lead to the recursion 
relations as 

dK^ dKz dgo , ^ 

^ = g^{3 -a-K-\l + A.(co + c,)/2K,)) - C,gl ^ = -C,gl8, (21) 

with g^i = Co 1(7, Cq 1 2 are positive constants and K^^ = ~ T. The recursion 

relations recover previous results on the Bose-glass and Vortex-Glass transitions, namely, 
for a = we have the Bose glass transition temperature = 3 without point disorder 

T8| . This critical value is decreased by both point disorder as found in Ref. 0, and by 



the exponent a. The Vortex-Glass transition KyQ = 2 is unmodified by columnar disorder 
or the exponent of the decay of the disorder correlations. Let us denote the dimensionles 
measure of disorder hj D: D = Ax{cq + ci)/2Kx- Now it was argued in Ref. p that when 
the point disorder is strong enough: D > 1, then the point disorder, measured by (70 > scales 
faster to strong coupling than the columnar disorder ((71), indicating a transition to a vortex 
glass phase. However this conclusion has to be viewed with some caution. Strictly speaking 
when deriving the scaling equations the coefficients (or equivalently the scaling dimensions) 
of both go and gi should have been small. Second, in typical real systems and numerical 
simulations the bare values of the forward and backward scattering are expected to be 
comparable. While there is no restriction on D in the theory, go has to be small, restricting 
the applicability of the results for D ^ 1 for typical systems. Therefore conclusions about 
a transition around D ~ 1 have to be viewed with some caution. 

Both of these problems are remedied by the present method. First we see that for a ~ 1 
the transition to the vortex glass phase happens at small disorder: D = 1 — a ^ 1. Second, 
in this region both scaling dimensions are indeed small, justifying a proper expansion. This 
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formulation thus constitutes a technically sound description of the transition to the vortex 
glass phase. 

However let us observe that the disorder exponent a scales in second order: for a < 
1 {6 > 0), a scales to zero, while for a > 1 (5 < 0) a scales to large values. Thus the system 
which started out on the Vortex Glass side of the transition 1>Z}>1 — a<^l will scale 
into the Bose Glass phase for a < 1. 

This flow into the Bose Glass phase for all a < 1 indicates an instability of the Vortex 
Glass phase towards the addition of arbitrarily weak columnar pinning. While the "weak- 
ness" is characterized by the critical exponent a being arbitrary close to the point-disorder 
limit 1, we expect that this result also translates to a critical value of c/p = 0. This result is 
also plausible as the same instability has been argued for single flux lines f^, and adding in- 
teractions will further localize the vortices. However, the Bose Glass phase itself is unstable 
to point disorder, although admittedly only on astronomically long length scales 0. The 
most likely picture emerging from the RG is therefore that the trajectories flow towards a 
strong coupling fixed point (s), which is characterised by simultaneous finite values of c and 
p. Because of the presence of columnar disorder, this feature would most likely translate to 
a finite localization length. 



B. Numerical Approach 

To test this scenario we have employed our Monte Carlo algorithm to investigate the 
phase boundary between the insulating and superconducting phases of the boson Hamilto- 
nian, corresponding to the Bose (large V) or Anderson (small V) glass and either the Line 
liquid or the Vortex Glass phase, respectively. No attempt was made to distinguish between 
the Vortex glass and the Line liquid phases since it has been difficult to find the Vortex glass 
phase from static properties [0. Rather, we instead measure the superfiuid density ps as a 
function of both boson repulsion V/t = 2(7ri^)^ and c/p to map out the Bose or Anderson 
glass phase boundary. 

Results of our runs are shown in Fig. 5 for the case of a 64 site chain with 40 bosons at a 
temperature (3 = 16. Values of ps are obtained by averaging over several realizations of the 
columnar disorder (fixed at c=2) and point disorder, although each run showed only slight 
variations in the superfiuid density and appears to be self averaging. Typical runs were over 
50,000 warm-up and 100,000 measurement sweeps. The bottom curve in Fig. 5 represents 
new measurements at lower temperatures of our previous results for ps |T^. The transition 
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from the superfluid phase to the locahzed one at large repulsions occurs at a slightly greater 
value than the Giamarchi-Shultz fixed point but is reduced compared to the runs which 
were performed at higher temperatures |T^. This may indicate that the localization length 
is nearing the chain size at this value of disorder. Further the localized phase at small 
repulsion (Anderson glass) becomes more pronounced at these lower temperatures. 

Turning on the point disorder shows that the superfiuid phase becomes enhanced, which 
is shown in the upper curve in Fig. 5 for c/p = 1. The Bose-glass transition now occurs at 
a larger boson repulsion, which is in agreement with the RG results from Eq. (21) and Ref. 

Further, the Anderson glass is even more adversely affected than the Bose glass phase 
as there is no evidence of a localized phase at small boson repulsion. It may be possible that 
the localization length is beyond the length of the chain at this temperature although runs 
performed at larger /3 = 32 do not support this. This would indicate that the Anderson glass 
is destroyed due to the de-phasing effect of the presence of random site energies changing in 



time, and is in agreement with previous results concerning Mooij anomalies |]T5| . In addition, 
we have performed runs at c/p = 0.1 and have found no evidence for a localized phase for 
any value of boson repulsion. 



V. DISCUSSION AND CONCLUSIONS 



In this paper we have presented the analysis of independent and interacting flux lines in 
the presence of competing point and columnar disorder. Our Quantum Monte Carlo simula- 
tions have verifled that for point disorder only the exponent C, characterizing line wandering 
changes from 1/2 to 2/3, in agreement with theoretical predictions. In the presence of 
columnar disorder a single flux line becomes localized. The localization length grows with 
the introduction of point disorder. It was argued in Ref. |^ that the localization length re- 
mains flnite but very large for large values of c/p. Our results are not inconsistent with this 
picture. We found clear evidence for a flnite localization length for c/p < 0.8, but for larger 
values of c/p the localization length crosses over to very large values, which are inaccesible 
with our method. 

For interacting flux lines we studied the same issues with the help of Renormalization 
Group and Monte Carlo techniques. We developed a variable exponent description of the 
disorder, which allows for a technically sound scaling treatment of the problem. The scaling 
trajectories indicate the instability of the Vortex glass phase for arbitrarily small values 
of columnar disorder towards a Bose glass, which is presumably characterized by the lo- 
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calization of the vortices. Furthermore as the single flux hne appears to be locahzed for 
arbitrarily small columnar disorder as well, it is likely that adding interactions between the 
hues will further enhance localization. Our Monte Carlo technique found a Bose Glass phase 
at medium values of the columnar disorder. The lack of flnding this localized phase extend- 
ing to very small values of c might be due to system size restrictions and the importance 
of rare regions. The localization length could thus be large but flnite. It was also shown 
that the Bose glass phase is unstable towards the addition of point disorder, but only on 
astronomically long length scales. Thus presumably the ultimate asymptotic behavior is 
governed by a fixed point with both c and p finite. As c > 0, the main physical characteristic 
of this phase is expected to be the localization of the vortex lines. 

We continue by noting that previously it has been implicitly assumed that the relevant 
quantity which characterises the disordered environment of flux lines is the ratio of the 
columnar to point disorder strengths, c/p. However, the recursion relations of Eq. (21) 
seem to prevent a reduction in the scaling to a single dimensionless coupling due to the 
presence of the second term in the recursion relation for gi, which predicts the reduction of 
the Bose glass transition temperature. Further, since the RG of Eq. (21) is limited to small 
values of the couphng constants, the second order terms could also lead to the conclusion 
that both c and p {gi and go) are relevant parameters, leading to a three-dimensional phase 
diagram in K~^,gQ and gi. This can be tested numerically and is a subject for future 
consideration. Lastly, further numerical simulations are necessary to to resolve the nature 
of the strong coupling fixed point (s). 

In conclusion, we have studied the behavior of free and interacting fluxlines in the simul- 
taneous presence of point and columnar disorder. Renormalization Group analysis indicates 
in both cases that the addition of even weak columnar defects localizes the flux lines. How- 
ever the physics seems to be described by a strong coupling flxed point, the nature of which 
has to be clarifled by numerical methods. To this end Quantum Monte Carlo simulations 
were also performed, which show such a localization behavior for medium and strong colum- 
nar disorder only. The possible reason for that may be that the localization length exceeds 
the system size in the case of weak columnar pinning. It was also demonstrated that the 
transition temperature to the Bose Glass is reduced by the presence of point disorder. Fur- 
ther work is necessary (for instance from dynamical quantities) to distinguish between the 
Vortex and line liquid phases. 
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FIGURES 

FIG. 1. The transverse wandering correlation function for different positions in imaginary time 
P along the boson world line for the case of no disorder, Eq. (14). A turnover from linear diffusive 
behavior occurs when the flux line begins to see the finite lattice size. 
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FIG. 2. The transverse wandering correlation function for different positions in imaginary time 
P along the boson world line for the case of two point disorder strengths. The slope of the line is 
1.33 ± 0.03. 
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FIG. 3. Transverse wandering correlation function at imaginary time slice /3/2 for p = 10, c 
(crosses); p = c = (diamonds); p = 10, c = 5 (triangles); p = c = 10 (x 's); p = 0, c = 10 (squa 
for a 32 site chain. 
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FIG. 4. Square of the inverse localization length as a function of p/c for a 32 site chain. 
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FIG. 5. Superfluid density ps for various repulsion strengths V/t for the case c = 2,p = 
(bottom hne), and c = p = 2 (top hne). 
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